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We investigate how a horizontal plane Poiseuille shear flow changes lin- 
t> ■ ear convection properties in binary fluid layers heated from below. The full 



linear field equations are solved with a shooting method for realistic top and 
bottom boundary conditions. Relevant characteristic growth exponents and 



the spatial structure of their associated eigenfunctions are evaluated for dif- 
^ ■ ferent perturbations of the conductive state. Through-flow induced changes 



of the bifurcation thresholds (stability boundaries) for different types of con- 
vective solutions are determined in the control parameter space spanned by 



^ . Rayleigh number, Soret coupling (positive as well as negative), and through- 

flow Reynolds number. We elucidate the through-flow induced lifting of the 
Hopf symmetry degeneracy of left and right traveling waves in mixtures with 
negative Soret coupling. Finally we determine with a saddle point analysis 
of the complex dispersion relation of the field equations over the complex 
wave number plane the borders between absolute and convective instabili- 
ties for different types of perturbations in comparison with the appropriate 

Ginzburg-Landau amplitude equation approximation. 
PACS:47.20.-k,47.20.Bp, 47.15.-x,47.54.+r 
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I. INTRODUCTION 

An externally imposed flow can influence the spatiotemporal behavior of dissipative 
structures growing in forced nonequilibrium systems. Examples are chemical and reaction- 
diffusion systems, biological problems, and the large variety of different hydrodynamic in- 
stabilities leading to pattern formation M. Here we theoretically investigate with a linear 
analysis of the relevant field equations the spatiotemporal properties of convection solutions 
that bifurcate out of the homogeneous, quiescent conductive state in a binary fluid layer 
heated from below. 

A lot of experimental @-[K|, analytical [pT-171, and numerical [18[-|211 activities have 



been devoted recently to investigate these primary convection patterns in the absence of 
through-flow ||22|| . They revealed a variety of bifurcation properties and spatiotemporal 
behavior that is much richer than that of stationary mirror symmetric roll patterns growing 
in the supercritical bifurcation of the standard Rayleigh Benard setup with a pure fluid, 
say water. In binary mixtures like, e.g., ethanol- water there occur stationary square and 
roll patterns but also symmetry degenerate left or right traveling convection waves. The 
bifurcation of the traveling wave (TW) solutions and of the stationary roll patterns can 
be super-, tri-, or subcritical relative to the critical heating rate. These different solution 
properties are controlled by the combination of thermal forcing — i.e. the Rayleigh number 
- and the strength ip [jl]] of the Soret coupling between temperature and concentration 
field. It is the concentration field that causes the rich structure formation behavior via its 
contribution to the buoyancy force field which drives convection. 

Without lateral through-flow several experimental and theoretical papers have addressed 
the linear convection properties |25|-p7|. The most comprehensive and most accurate results 



were obtained in the more recent numerical work [pT| -|37[|. In this work we determine how a 
lateral through- flow changes the structure and dynamics of convection fields at onset. We 
evaluate relevant characteristic exponents and the spatial structure of their associated eigen- 
functions for different perturbations of the conductive basic state and determine the through- 
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flow induced changes of the bifurcation thresholds for different types of convective solutions. 
These thresholds are defined by vanishing real parts of the characteristic exponents. They 
mark stability boundaries of the conductive state against different convective perturbations 
in the control parameter space spanned by the Rayleigh number Ra, the Soret coupling ijj, 
and the through-flow Reynolds number Re. We have numerically solved the full linear field 
equations subject to realistic boundary conditions using a shooting method for perturba- 
tions with wavevectors parallel to the through-flow. The solutions for other wavevectors can 
be obtained from the former by a straightforward symmetry transformation. Results were 
obtained for positive as well as for negative ip. For ip < the most important result is that 
the through-flow lifts the Hopf symmetry degeneracy of left and right traveling waves at 
Re = 0: frequencies, bifurcation thresholds, and structural properties of the two waves are 
changed dramatically. Mixtures with more negative ip require a larger Re for the changes 
to reach a comparable relative size. For sufficiently large Re the lowest relevant bifurcation 
threshold of binary mixtures with any ip asymptotically approaches the critical Rayleigh 
number Ra c (Re,ip = 0) of a pure fluid with imposed through-flow. Then the externally 
imposed shear flow eliminates the Soret-induced coupling effects between the convective 
concentration field and the other fields by suppressing vertical convective transport of Soret 
driven concentration perturbations. 

Our paper is organized as follows. In Sec. |J we describe the system. In Sec. |T| we review 
the linearized equations for perturbations of the conductive state, their boundary conditions, 
the eigenvalue problem for the characteristic exponents, relevant symmetry properties, and 



the behavior for small through-flow rates. Sec. IV contains our results concerning bifurcation 



properties — stability thresholds, wave numbers, frequencies, and eigenf unctions - - for 
negative and positive Soret coupling ip as functions of the through-flow Reynolds number 
Re. In Sec. |V| we compare borderlines between absolute and convective instabilities obtained 
for different types of perturbations from the field equations with results from the Ginzburg- 
Landau amplitude equation. The last section gives a brief summary of our work. Appendix ^ 
contains details of our shooting method. There we also describe our procedure to find saddle 
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points of the dispersion relation of the field equations in the complex wave number plane. 
Appendix || presents results obtained from a variational calculus. 



II. SYSTEM 



We consider a horizontal layer of height d of a binary fluid mixture in the homogeneous 
gravitational field g = — ge z that is directed downwards. A positive temperature difference 
AT is imposed between the lower and upper confining boundaries, e.g., via highly conducting 
plates in experiments. The associated Rayleigh number is 



where k is the thermal diffusivity and v the kinematic viscosity. The thermal expansion 
coefficient a and the solutal expansion coefficient f3 follow from a linear isobaric equation of 
state for the total mass density 



for small deviations of the temperature T from its mean T and small deviations of the 
solute's mass concentration C from its mean Co- 

An externally applied lateral pressure gradient drives a through-flow in x direction. The 
resulting mean lateral flow velocity U determines the through-flow Reynolds number 



We investigate here the parameter regime < Re < 1. With d ~ 0.5 cm, v ~ 0.01 cm 2 /s 
(H 2 0) the maximal averaged through-flow velocity is then U ~ 0.02 cm/s, i.e., 1.2 cm per 
minute. 




(2.2) 




(2.3) 



A. Equations 



To describe this system we use the balance equations for mass, momentum, heat, and 
concentration in Oberbeck-Boussinesq approximation |]23|j39| 
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V • u = (2.4a) 
(9( + u-V)u = ffV 2 u - Vp 

+ a(T-T + C-C )e z (2.4b) 

( + u • V ) T = V 2 T (2.4c) 

(3t + u- V)C = LV 2 (C-^T) . (2.4d) 

Here u = it + v e y + w e z is the velocity field. We reduce lengths by d, times by d 2 / n, 
the effective pressure p by K 2 /d 2 , temperatures by K,u/(agd 3 ), and the concentration field 
by Kz//(/3g<i 3 ). Then the two material parameters Prandtl number 

a = - 2.5 

K 

and Lewis number 

L — — (2.6) 

K 

appear with D being the concentration diffusion constant. Furthermore, there enters the 
separation ratio 

* = (2-7) 

that measures the strength of the linear Soret coupling between concentration and temper- 
ature field via the thermodiffusivity kr- 

B. Conductive state 



For small Ra, Re a laterally homogeneous solution of (|2.4|) is stable that describes a 
conductive state without vertical convective flow. It is a combination of plane horizontal 
Poiseuille flow 



Ucond = U (z) e x = a Re P{z) e x 
P(z) = 6z(l-z) 



(2.8a) 
(2.8b) 
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and a diffusive temperature field 

T cond = T + Ra(^-z) (2.9) 
that enforces via the Soret effect a diffusive vertical concentration stratification 

C cond = C -Raij{ l --z) . (2.10) 

III. CONVECTIVE PERTURBATIONS 

Here we briefly review the linearized equations for perturbations of the conductive state, 
their boundary conditions, the eigenvalue problem for the characteristic exponents, relevant 
symmetry properties, and the behavior for small through-flow rates. 

A. Linearization around the conductive state 

1. Equations 

The basis for our linear analysis of convective perturbations of the conductive state 
described in Sec. [II B| are the linearized field equations 



( UV 2 - d 2 z U )d x w = a (81 + d 2 y ) {9 + c) (3.1a) 

(d t - V 2 + Ud x )6 = Raw (3.1b) 

(d t - L V 2 + Ud x )c = ipRaw - Li/jV 2 6 (3.1c) 

for the deviations 

= T — T con d, c = C — C conc i (3-2) 



from the conductive state (|2.9| , [2.10| ). Here w is the vertical velocity field that vanishes in 



the conductive state. To derive ( |3.1a| ) we have applied twice the curl operator to equation 



( [2.4b|) using ( |2.4a| ). Note that the Poiseuille flow profile U(z) of the conductive state enters 



into Eqs. ( |3.1|) making them non-autonomous. 
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2. NSI boundary conditions 

We consider the horizontal boundaries to be perfectly heat conducting and rigid with no 
slip and vanishing vertical concentration transport. These, so called, NSI (no-slip, imper- 
meable) conditions impose 

= w = d z w = d z ( = at z = 0, 1 . (3.3) 

Here we have introduced the combined field 

C = c-# . (3.4) 

Since the concentration current at the no-slip boundaries is purely diffusive the condition 
d z ( = ensures impermeability of the horizontal boundaries. Laterally we assume the 
system to be unbounded. 

3. FS boundary conditions 

As an illustrative special case let us consider for the moment free slip (FS) horizontal 
boundaries with a shear-free plug flow profile, Ups{ z ) — U. In this idealized situation the 
effect of through-flow in Eqs. (|3.1|) can be transformed away by a Galilei transformation to 
a system that comoves with the vertically constant plug flow velocity U = aRe. Thus, the 
stability properties of the FS conductive state are not changed by the horizontal plug flow. 
The stationary and oscillatory marginal stability curves of the mixture remain the same, 
only the characteristic exponents acquire an additional imaginary part of size ik x U. This 
holds for FS horizontal boundaries irrespective whether they are permeable or impermeable 
to the concentration field c. 

B. Eigenvalue problem 

The general solution of the perturbation equations ( |3.1| ) can be written as a superposition 
of plane-wave perturbations with lateral wavevector 
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(3.5) 



The plane-wave solution ansatz for the fields 



(3.6) 



reads 



(3.7) 



with a complex characteristic exponent 



s = 9fts + i = 7 — iuj 



(3.8) 



and complex ^-dependent amplitude functions $ = (w,9, Q. Inserting the ansatz (|3. 7|) into 
the field equations (|3.1| ) yields the 3x3 linear eigenvalue problem 



{C + sM) ${z) =0 



(3.9a) 



for the eigenvalues s and eigenvectors $ with 



C = £ (0) + iak x Re£ {1) 



(3.9b) 



Into 



( -a(d 2 -k 2 ) 2 a(l + ^)k 2 ok 2 ^ 



£(0) 



-Ra 




£(1) 



A; 2 - <9 2 
</>(<9 2 -£; 2 ) L{k 2 -d 2 z ) 



P (d 2 - k 2 ) - d 2 P 
P 

OP 



(3.9c) 



(3.9d) 



enters the vertical profile P(z) ( |2.8b ) of the Poiseuille through-flow and its second derivative. 
Finally 
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d 2 - k 2 
1 




(3.9e) 



In the absence of through- flow C reduces to C^°\ 

Due to the boundary conditions the eigenvalue spectrum is discrete. We are interested 
in the three characteristic exponents Sj (j = 1, 2, 3) whose growth rates jj are closest to 
zero and whose eigenfunctions &j(z) have no nodes other than at the horizontal boundaries 

2 = 0,1. 

C. Symmetries 

The solution of ( |3.9|) , i.e., eigenvalue s and eigenfunction $ depend on the material 
parameters a and L, the control parameters ip, Ra, Re, and on the lateral wavevector 
k. Note first of all that the dynamics of perturbations with wavevectors perpendicular to 
the through-flow are not changed by the latter, since for them the contribution from 
vanishes when k x — 0. 

1. Squire transformation 



Since k and Re enter into (|3.9| ) only as k 2 and k x Re the dependence of the functions 
f — s, $ on k and Re is 



/ = / (k 2 , k x Re 



(3.10) 



Using this behavior the Squire transformation 



f(k 2 x + k 2 y} k x Re) =f(kl,k x Re) 
kl = k 2 x + k 2 y ■ Re= k ^Re 

krr 



(3.11a) 
(3.11b) 



relates the functions / for a wavevector with arbitrary components k x , k y to the functions 
that have been determined for wavevectors k = k x e x in through-flow direction and Reynolds 
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numbers Re. Therefore, we shall consider in the remainder of this work only k y = pertur- 
bations with wavevectors k = k x e x that are parallel or antiparallel to the through-flow. For 
Re = the Squire relations ( |3.11|) reflect the horizontal rotational symmetry of the system 



in the absence of through-flow. 

2. Reverting the through-flow direction 

Upon reverting the flow direction, i.e., under the operation Re — > —Re the set {sj, $.,} 
of eigenvalues and eigenfunctions transforms into each other since the balance equations 
( |2.4|) are invariant under the parity operation (x,u) — > —(x,u) with u being the velocity 
field in x-direction. The transformation behavior of {sj, $.,•} follows explicitly from the fact 
that the linear operator L entering the eigenvalue equation (|3.9|) transforms as 



= £*(<lx) ; <lx = k x Re (3.12) 

under Re — > —Re with the star denoting complex conjugation. Here we do not display the 
other arguments of £ that remain unchanged. We use ( |3.12|) in the complex conjugate of 
Eq. 



[£(-?») + s*(q x )M] = . (3.13) 

Thus, if Sj(q x ) with &j(q x ) are solutions of (|3.9|) so are 

Sj(q x ) = s*(-q x ) with S^fc) = (3.14) 

Now, the nondegeneracy of the eigenvalue problem implies that the two sets {sj,$j} and 
{sj, $j} are the same. We find that one eigenvalue, say j — 3, does not change — s 3 = s 3 
and $ 3 = $ 3 — so that according ( 3.14j) 



73(-?<c) = 73(?<n) (3.15a) 
-U3{-q x ) = uz(q x ) (3.15b) 
= $ 3 (q x ) . (3.15c) 
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The other two eigenvalues and eigenf unctions are crossrelated to each other — s% — S2, $i 
$2 and §2 = Si, $2 = $i — so that according to fl3.14| ) 



7i(-&) = 72(g*) 

-^i(-gx) = ^2(9*) 



(3.16a) 
(3.16b) 
(3.16c) 



For i?e = the relations (|3.15| , |3.16| ) reflect the facts that in the absence of through- 
flow (i) the solutions depend on k 2 x only, (ii) one eigenvalue is real, say u 3 = 0, with real 
eigenfunction $3 = $3 that describes nonoscillatory dynamics, and (iii) the other two are 



a complex conjugate pair 



7i 



72, ^1 



-u 2 with $1 = $ 5 



that describes the 



growth/decay of symmetry degenerate left and right traveling waves. This symmetry is 
broken by the through-flow. If for Re = all three eigenvalues and eigenf unctions are real 



then they fulfill for finite Re relations like ( |3.15| ). 

Since ( |3.15| , |3.16[ ) relate the solutions for negative k x or Re to those with positive k x , Re 
it suffices to investigate the three eigenvalues Sj and eigenfunctions $j for positive k x , Re 
only. 



D. Expansion for small Reynolds numbers 

It is instructive to see how the Re = solutions of the eigenvalue problem 

( p.9|) evolve upon switching on the through-flow. The qualitative behavior can be studied 
analytically for small Re via an expansion in the parameter 

r] = aq x = ak x Re (3.17) 

that appears explicitly in the linear operator £ ( |3.9b|) 

C = C {0) +ir]C {1) . (3.18a) 

So we expand 
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s . = sf ) +rjsf ) +0(rf i ) 
$ i = $f +77$f+0(?7 2 ) . 

Inserting Q3.18|) into (|3.9a|) yields in order 77 the equation 

'(0) , JO) A A \ <i,"> 



(3.18b) 
(3.18c) 

(3.19) 



which is solvable under the condition 



(o)t 



;(o) 





Here $^°' ) * is the solution of the adjoint equation 

of the zeroth-order eigenvalue problem 

(£(°) + sf JW) $f = 



(3.20) 



(3.21) 



(3.22) 



s,° is the complex conjugate of the original 



for the eigenvalue Sj. Note that s^- 
eigenvalue s^. The Fredholm alternative ( |3. 20| ) leads to the first-order correction for Sj 



where 



Pj 



4" = - 













(3.23) 



(3.24) 



is determined by normalized "matrix elements" of through-flow perturbation "operators" 
containing the profile P(z) between the zeroth-order eigenf unctions = (wf ,0j ,Cj 
Thus, the small- Re expansion yields the following results for the eigenvalues Sj = 7^ 



jj = 7,- 0) + ?7 ^Pj + 0(r) 2 



(3.25a) 
(3.25b) 



Uj = uf +r ] Mp J +0(r ] 2 ) . 

For the subsequent discussion we use the fact that the "operators" entering into the 
"matrix elements" of Eq. (|3.24|) are real. 
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1. Stationary perturbations 

Let us call a perturbation stationary, for shorthand, if the characteristic exponent sf^ 
in the absence of through-flow is real. Then the corresponding eigenfunction is real as well, 
= which implies that also pj ( |3.24| ) is real. Thus, with Qpj = 0, one obtains from 

7j . = 7 j 0) + 0(i?e 2 ) (3.26a) 
UJj = ujf + ak x Re % + 0(Re 2 ) . (3.26b) 

The frequency grows linearly with k x Re while the growth rate is an even function of k x Re for 
these stationary perturbations as shown schematically in Fig. |1] for the eigenvalue labelled 
by j = 3. 



2. Oscillatory perturbations 

We call the perturbations oscillatory that are described by the two eigenvalues (say, j = 
1,2) which in the absence of through-flow form a complex conjugate pair, = s^*. The 
corresponding eigenfunctions are complex conjugates of each other, $2°'' = ^f^*, according 
to ( |3.16| ) which implies P2 = p*. With 9ffy>2 = 3^Pi and Qp2 = —^Pi one obtains from ( |3.25| ) 
the following relations 

71 = 7j 0) ± ak x Re Spi + 0(Re 2 ) (3.27a) 
ui = ±c;i 0) + a^ReUpx + 0(Re 2 ) . (3.27b) 

Since with the plane- wave perturbation ansatz ( |3.7|) the TW phase velocity u/k x should be 
increased for positive k x , Re by the thro ugh- flow one can expect without having performed 
an explicit calculation that 3lpi in ( |3.27b| ) should be positive — cf. Fig. [I]. And consequently 
both frequencies id± \ are shifted upwards by the same amount ak x Re^tp\. On the other 
hand, to decide whether the growth rate of the left or of the right traveling wave is increased 
or decreased requires an explicit calculation of In any case, however, Eq. (|3.27a|) 
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predicts that for small through- flow the symmetry degeneracy, = 72 , of the growth rates 
is lifted by a symmetric splitting that increases linearly with Re as indicated schematically 
in Fig. ^ for 71 and 72. 

For idealized FS boundary conditions the frequencies ujj = + a k x Re behave as in 
Fig. [1]. However, the growth rates are independent of Re so that the degeneracy 71 = 72 
is not lifted by the FS plug flow. 

IV. STABILITY AND BIFURCATION PROPERTIES 

Here were present for negative as well as for positive Soret coupling if) the Re- dependence 
of critical properties: stability thresholds, wave numbers, frequencies, and eigenfunctions. 
These results have been obtained numerically by a variant [[|!| (cf. Appendix Kl for a short 



description) of a standard shooting method that has previously been used |30| to determine 
stability properties of binary mixtures in the absence of through-flow ||36|| . In order to check 
these results — in particular some of the unexpectedly strong and peculiar changes with 
Re — by an independent method we have performed a variational calculation described in 
Appendix [TJ. 

A. Notation 

We introduce Rayleigh numbers and wave numbers 

(4.1a) 



Ra c (Re = Q,if) = 0) 
k 

k= k c (Re = 0^ = 0) (4 ' lb) 
that are reduced by the critical ones of a pure fluid (if> = 0) in the absence of flow (Re = 0) 

Ra c (Re = 0, if) = 0) = 1707.76 (4.2a) 
k c (Re = 0,ip = 0) = 3.11632 . (4.2b) 
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As explained in Sec. [Ill C| , it suffices to consider wave vectors k = k x e x in through-flow 
direction with positive component k x = |k| = k. The spatiotemporal behavior of pertur- 
bations with other vectors follows with the symmetries of Sec. [Ill C[ We determine the 
evolution of the three relevant eigenvalues Sj — jj — i uij (cf. Sec. [Ill B| ) upon increasing Re 
from Re = up to about 1. In particular we evaluate the critical parameter combinations 
for which each of the three growth rates jj first passes through zero when increasing r. 

The so obtained critical quantities are marked by a subscript c and in addition by a su- 
perscript S, U, D that replaces the running index j of the three eigenvalues used in Sec. fH\ . 
The superscripts S, U, and D identify the critical perturbation behavior, e n k cx-u} c t) ^ j n ^ e 
limit Re —>■ 0. Eigenvalues for which uj c {Re — > 0) = are marked by 5* since these pertur- 
bations are stationary for Re = 0. Eigenvalues for which u c (Re — > 0) is positive (negative) 
carry the superscript D (U) since they characterize for Re — *> perturbations which prop- 
agate in downstream (upstream) direction. We should like to stress again that the cases 
5* ("stationary"), D ("downstream"), and U ("upstream") characterize the perturbations 
in the limit Re ^0 — see also Fig. |1} In general all critical frequencies are finite in the 
presence of through-flow. However, for a special value of Re one has uj^ = while cuf and 
ujD are positive (cf. Sec. |1V lj| and Sec. [IV C]). 



B. Effect of through- flow on the oscillatory instability 

We present in this subsection critical properties of a binary mixture like water-ethanol 
with L = 0.01, a = 10, and separation ratio ip = —0.1 as a representative case for a moder- 
ately negative Soret coupling. For these parameters the nonlinear solution of stationary con- 
vection in the absence of through-flow, Re = 0, is already disconnected from the conductive 
state since the stationary bifurcation threshold has moved already at ip^ = —L/(l + L) 
to infinity PTfl . However, at = = 1.1200 there is for Re = an oscillatory threshold 
into symmetry degenerated left and right — or in our language upstream and downstream 
- propagating traveling waves with critical wave numbers k® = k^ = 1.0022 and critical 
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Hopf frequency = —u^ = 6.4659. 

In Fig. we show the variation of these two critical thresholds, wave numbers, and fre- 
quencies with increasing through-flow Reynolds numbers. In each case the full (dashed) 
line represents the upstream (downstream) critical quantity. Initially, for small Re, the 
bifurcation threshold r^(r^) is depressed (enhanced) by the through-flow so that the con- 
ductive state is destabilized (stabilized) against convective perturbations traveling upstream 
(downstream) in comparison to the symmetry degenerate Hopf bifurcation threshold with- 
out through-flow — see also the schematic variation of the growth rates in Fig. [I]. So at 
very small Re it is upstream traveling wave convection that grows first when increasing the 
Rayleigh number quasistatically. 

However, while increases monotonically with Re — first linearly and then quadrat- 
ically - - the initial linear downwards shift of changes at larger Re to a precipitous 
increase and a subsequent flattening. Thus, in the shown .Re-range the two critical curves 

and have two intersections giving rise to bistable behavior of perturbations there. 
Note, however, that the bistable upstream and downstream TW perturbations can not be 
superimposed linearly to a standing wave since their wave numbers k]! ^ k® differ and 
furthermore ^ So in the .Re-interval between the bistable intersections of r]! and 

r,? downstream propagating convection waves grow first while outside this interval at small 
Re and large Re upstream TW convection bifurcates first out of the conductive state. 

Considering the critical wave numbers k^ and of Fig. [|b it should be noted that their 
variation is very small — less than 2% — and that /c^ 7 < k^ in the .Re-range of Fig. |2|. This 
behavior holds also for the other Soret coupling strengths ip = —0.25, —0.01, and — 0.001 
that we have investigated [fl3[ . 



Somewhat unexpected to us is the non monotonical variation of k^ and also of dur^ / dRe 
(Fig. 0c) in the interval below Re = 0.5 where uj^ changes sign and where r^ 7 shows its 
strong increase. To check that this variation is not a numerical artifact of our shooting 
algorithm we have performed a stability analysis with a variational approximation being a 
fundamentally different method. The variational results presented in Appendix [B] also show 
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the peculiar variation of with Re obtained from the shooting method thus supporting 
the latter behavior. 

The critical frequencies uj^ and shown in Fig. |2|c are practically linear functions of 
Re and in this respect similar to the frequencies of the idealized FSP system. They start 
at zero through-flow with the Hopf values l^ -*] and — |c<4°'|, respectively, and they can be 
very well approximated by the first-order result ( |3.27b|) of the low- .Re expansion 

Uc = o4 0) + a *4 0) Kpi Re (4.3) 

with ak^dtpi ~ 41.9 for a = 10. Comparing this rate of change du c /dRe ~ 41.9 with 
results for other separation ratios including the pure fluid case E3] one finds only very small 



deviations JJ3J. Obviously dtpi ( |3.24| ) depends only weakly on the Soret coupling ip. 

Note that for Re > |o;^|/41.9 both critical frequencies, and , are positive. Then 
the phase velocities, v c = w c /k c , of the two different critical TW's are in positive x-direction 
in the laboratory system, i.e., in through-flow direction. However, is always smaller - 
by about 2|v^| — than v®- Only for Re < |u;^|/41.9 is the phase velocity negative, i.e., 
opposite to the through-flow. So the wording "upstream propagating perturbations" that 
we are using in this work does not necessarily imply that the phase velocity of such a TW is 
negative in the laboratory frame. It would be negative in a frame that is moving in through- 
flow direction with a conveniently defined mean lateral velocity like, e.g., v = \{v® + v^). 

C. Bifurcation thresholds at negative ip 

In Fig. |] we show the bifurcation thresholds (full lines), (dashed lines), and r~ 
(dotted lines) as functions of Re for a few characteristic negative Soret couplings if). 

1. r u c (Re^) 

The typical shape of the stability curve r^ 7 that is displayed in Fig. ^]a for ip = —0.1 does 
not change much for other negative separation ratios: As a function of Re (full lines in 
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Fig. ^|) decreases for small Re, develops a minimum where u>^ goes through zero, steeply 
increases thereafter, and finally flattens asymptotically towards r^(Re,ip = 0) at large Re. 
Thus a sufficiently large through-flow eliminates the Soret induced coupling effects between 
concentration field on one side and temperature and velocity field on the other side: The 
bifurcation threshold r^(Re,ip) approaches for any Soret coupling ijj the pure fluid stability 
boundary rf (Re, ip = 0) at large Re. For small if), e.g., at ip = —0.001, the stability boundary 
lies always below r;? while for larger (see, e.g., ip = —0.01) there are two intersections 
of the curves r,F and with r? < rY in between — cf. the related discussion in Sec. ffVB 



2. r^(Re,ip) 

The bifurcation threshold (dashed lines in Fig. ^) always increases monotonically with 
the through- flow strength. The initial slope drf? jdRe increases somewhat with decreasing 
|^| . For ip = —0.001 and ip = —0.01 the stability curves and collide in the .Re-range 
displayed in Fig. [3| This property is elucidated in paragraph 4 further below. 

3. rf(Re,ip) 

In pure fluids, ip = 0, the bifurcation threshold r^(Re,ip = 0) (lowest dotted curve 
in Fig. ^|) slightly increases with growing Re []42 |. In binary mixtures with negative Soret 



coupling, on the other hand, the stationary threshold rf gets very strongly depressed by a 
small through-flow — see the dotted curve for ip = —0.001 that starts at rf(i?e = 0,ip = 
-0.001) = 1.1816. 

In the absence of through-flow, Re = 0, the threshold rf rapidly increases with \ip\ and 
diverges at ip^ = —Lj (l + L) = 0.0099 for L = 0.01. Beyond this Soret coupling the solution 
branch of stationary nonlinear convection is disconnected from the ground state solution as 
r^(Re = 0, ip < ip^) = oo. A small but finite through-flow, however, moves the threshold 
rf down to finite values: The dotted curve for in the inset of Fig. [5] for ip = —0.01 < i/j^ 
shows (i) that rf = oo below a finite Re^ ~ 0.019, (ii) that rf is finite for Re > Re^, 
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and (in) that steeply drops down for Re > Re^. The Reynolds number Re^ where 
diverges grows with increasing \%p \ — a stronger Soret coupling requires a larger through-flow 
to move the bifurcation threshold rf from infinity to a finite value. 

4- Collision of the and rf stability boundaries 

With increasing Re the bifurcation threshold and approach each other. The former 
decreases rapidly and the latter increases with Re and they almost coalesce in the Re — r 
plane of Fig. This behavior is most easily understood by investigating how the relevant 
eigenvalues s = 7 — iu vary with r and Re. To that end we show in Fig. |] 7(7") and u(r) of 
the two relevant eigenvalues for a representative Soret coupling ip = —0.01 at Re = (thick 
curves) and at Re = 0.4 (thin curves) for a fixed wave number k = 1. The real part of the 
third eigenvalue is always negative in the parameter range of Fig. |] and thus irrelevant for 
the following. Fig. |5| shows in a schematic way the motion of these two eigenvalues in the 
complex s-plane with increasing r for Re = (thick curves) and for a small Re 7^ (thin 
curves) . 

Let us consider first Re = 0. Then there is for small r a complex conjugate pair of 
eigenvalues (7^ = j d ,uj u = — u D ) that produce the symmetry degenerate Hopf bifurcation 
(at r D = r u « 1.021 in Fig. f|b) when 7^ and 7 D pass simultaneously through zero. This 
situation is marked in Fig. |] and Fig. |5| by thick upwards and downwards pointing triangles. 
With increasing r the real parts 7 grow and the frequencies u approach zero and the pair 
of eigenvalues meets in the complex s-plane of Fig. [5] on the real 7-axis (i.e. in Fig. |] at 
r m 1.231). Then, at larger r, this pair of real eigenvalues splits and moves apart along the 
real 7-axis (cf. thickly dotted SI and 5*2 curves in Fig. |] and [|). So we have a transformation 
of two oscillatory eigenvalues (U and D) into two stationary ones (5*1 and 5*2). For Re = 
the two thickly dotted branches 7 s1 and 7 52 in Fig. ||b remain above zero, i.e., SI does not 
reach the imaginary axis in Fig. |5|. 

Now for finite Re the symmetry degeneracy of the Hopf eigenvalue pair is lifted. The 
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r- value where 7^ goes through zero (A at r u « 1.05 in Fig. ^0) differs from the one where 
7 D = (v a t r D « 1.112 in Fig. |4]b) and the frequencies cj 17 and c<j d (thin lines in Fig. ^a) 
are shifted upwards. Thus the Re = pitchfork topology of the eigenvalue paths in Fig. |5] is 
perturbed. Moreover, the eigenvalue branches become disconnected (thin lines in Fig. |) as 
the thick eigenvalue branches of Fig. |] are deformed by the trough-flow into the thin ones - 
the arrows in Fig. f|b indicate the deformation directions. In particular the lower 7 51 branch 
of Fig. |]b and similarly the left moving SI branch in Fig. |5| goes at a sufficiently large Re 
through zero at the open circle (o) in Fig. f| and [| thereby producing an S'-instability for 
k — 1 at o, i.e., at a finite value of r s (« 1.296 in Fig. |J). Coming from 00 the intersection 
o in Fig. |]b has moved with increasing Re to the left to the finite value r s . By increasing 
Re further the thin D / Sl-curve in Fig. |]b is pushed downwards towards smaller 7, i.e., to 
the left in Fig. [5| Thereby the zero crossings at r D (xj) and r s (o) move together and vanish 
simultaneously at a particular .Re-value. Thereafter there is only the zero crossing of 7^ at 
r u (A). This is in principle what happens when the and rf curves in Fig. |3] at ip = —0.001 
and ijj = —0.01 approach each other with increasing Re. 

5. Opening of a wave number gap in the D — S marginal stability curves 

The above described merging of the zeros of the 7 13 / 51 curve that occurs in Fig. f|b for 
k — 1 slightly above Re = 0.4 corresponds to the opening up of a gap in the marginal 
stability curves against D and S perturbations. This scenario is documented in the k — r 
plane of Fig. |6|. There we show with gray-scales the height distribution of 

yvsi i n 

the k — r 

range where 7 > 0. Thick lines labelled by 7 = are marginal stability curves. In the white 
parts of Fig. |6| >y D / sl is negative and the dashed lines indicate 7-isolines for 7 < 0. The third 
eigenvalue j u / S2 being positive — cf . Figs. |] and [5] — is not shown. The open triangle and 
the circle in Fig. |]a mark the zeros of ^ D I S1 at Re = 0.4 that are shown in Fig. [|b by the 
same symbols. 

Upon increasing Re the eigenvalue >y D / sl decreases. Thus, the "mountain landscape" of 
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ryD/si somewhat globally "sinks" down. Thereby the 7 = isolines come together in the 
k — r plane of Fig. ^ and since the " mountain ridge" of ^ D / Sl does not have constant height 
(cf. gray-scales in Fig. ||a) the 7 = curves are connected into two tongues (Fig. |]b) that 
are separated by a wave number gap in which ^ D / S1 is negative. Increasing Re further the 
gap widens and the k — r regions with 7 > between the 7 = isolines narrow down as the 
latter move away from k = 1 towards larger r. 

We have difficulties to resolve this behavior of the 7 = isolines with our shooting 
method. We therefore show in Fig. |3] and later in Figs. |] and |1^ the minima and r~ of 
the marginal stability curves only before the wave number gap opens. The ending of r® and 
rf in Figs. [3], || and ^ should therefore not be interpreted as a termination of bifurcation 
branches: after the opening up of the gap the tongue shaped stability curves move towards 
larger r and with them their minima. 



D. Bifurcation properties at positive ip 

In the absence of through-flow there is only a stationary bifurcation threshold rf (Re = 
0, ip) at ip > that strongly drops from r^(Re = 0, ip = 0) = 1 towards zero when increasing 
ip. Switching on the through-flow has the overall effect of increasing rf towards rf (Re, ip = 0) 
as can be read off from Fig. |7|. Thus the lateral flow stabilizes the basic state by eliminating 
in the so called Soret regime |§ the convectively induced concentration homogenization. 
Note that already a very small through-flow has a dramatic stabilization effect: rf increases 
very strongly for small Re. Similarly the critical wave number fcf (Fig. 0b) approaches with 
increasing Re the pure-fluid value k^(Re,ip = 0). 



E. Bifurcation surfaces in r — Re — ip space 

To give an impression of the form of the three critical surfaces in the r—Re—ip space where 
U, D, and S convection patterns bifurcate out of the conductive state we combine in Fig. |8| 
in a 3-D plot the .Re-dependence of the bifurcation thresholds (thin full lines), r® (thin 
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dashed lines), and rf (thin dotted lines) presented so far together with their -^-dependence 
for Re = (thick lines). The bifurcation surfaces r^(Re,ip) and r^!(Re,ip) emanate for 
negative Soret coupling ip < out of the degenerate Hopf threshold line r®(Re = 0, ?/>) = 
r^(Re = 0, -0) (thick dashed and full line) and split apart when the through-flow is switched 
on. Upon increasing Re further r^ 7 gets indented slightly. On the other hand, curls up 
and comes very close to the surface r^(Re, ip) that is strongly bent down by the through-flow 
for ip < 0. See Sec. [IV C 5| for a discussion of the further fate of these bifurcation surfaces. 



The physically relevant, i.e. lowest lying, surface at larger Re is r^(Re,ip) (thin full 
lines). For any ip < it asymptotically approaches with increasing through- flow the ip = 
stability threshold, i.e., r^!(Re,ip) — > rf(i?e, ijj = 0). Thus a sufficiently large through-flow 
effectively eliminates the influence of any Soret coupling between convective concentration 
field and temperature and velocity fields (the diffusively induced concentration stratification 
in the conductive state ( [2.10|) , on the other hand, is not altered by the lateral shear flow 
). The bifurcation thresholds of the mixture approach that one, r^(Re,ip = 0), of the pure 
fluid. This also holds for mixtures with positive Soret coupling — cf. the thin dotted line 
for rS(Re,ip) at if) > 0. 



F. Structure of critical convective patterns 

Here we show how the spatial structure of the critical field deviations from the conductive 
state changes with increasing Re. To that end we have evaluated the complex eigenfunctions 
$(a;, z, t) = <£>(z) e l i k ^ x -^t) at ^he critical thresholds rf, r v c , and rf. Thus the critical 
convective fields w,6,c have the form 

w(x, z, t) = \w(z)\ cos [k c x — u c t — (p w (z)] (4.4) 

and similarly for 9 and c. Being solutions of complex linear equations we choose the abitrary 
complex scaling constant by fixing the modulus of the convective temperature field in the 
middle of the layer, 
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\e(z= 1 -)\ = l } (4.5) 
and by fixing the vertical mean of the phase f w (z) of the vertical velocity field to zero 



i 



o 



dz <f w (z) = 0. (4.6) 



We present gray-scale contour plots of the fields w, 8, and c in the vertical x — z cross 
section of the fluid layer with white (black) denoting large (small) values. Each of these 
figures H 11], [13| has nine contour lines denoting the fractions ±n/5 of the maximal field 



values with n = 0, 1, 2, 3, 4. In addition we present in figures 10|, [T|, [14] vertical profiles of 



the moduli and phases of the complex field amplitudes w,6, and c. 

1. Propagating patterns for Re = 
First we briefly recall the critical TW field structure (top part of Fig. in the absence 



of through-flow |3l||3"0| -|3"2"f for a relatively large Soret coupling ip = —0.25. The TWs of 
velocity, temperature, and concentration are vertically not plane but their phases show a 
vertical variation that is largest for the concentration wave (Figs. |9| and |i~0|) . The lateral 
location of a concentration surplus (white ellipse in the c-field of Fig. ^) phase-lags by about 
a quarter wavelength behind the lateral position of vertical downflow (black ellipse in the 
w-field). Thus, since the w-field advectively transports concentration surplus (deficiency) 
from the Soret induced alcohol rich top (poor bottom) boundary layer into the bulk fluid this 
feeding mechanism between boundary layer and bulk is laterally phase shifted by roughly 
A/4 in the propagating wave. Also the crest (valley) position of the temperature wave 
phase-lags — albeit by a smaller amount — behind the lateral location of maximal vertical 
upwards (downwards) flow which advectively feeds the bulk with warm (cold) fluid from the 
warm bottom (cold top) region. This phase-lag agrees quite well with the result 
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obtained |16| from a Galerkin model. 



As a result of the smallness of the Lewis number L, i.e., of the diffusive concentration 
transport the c-field shows characteristic boundary layer behavior near the plates where 
advection decreases to zero — see the variation of the modulus \c(z)\ and of the phase tp c {z) 
in Fig. [H]. Whithin the boundary layers \c(z)\ is suppressed and the phase-lag of the c-wave 
behind the w-wave is significantly enlarged. 

2. "Downstream" patterns 

Upon turning on the through-flow the phaselines of the downstream propagating patterns 
get bent further (cf. Fig. |9| and in particular the right column of Fig. [10]) and the phase 



differences between the different waves increase near the plates. The vertical moduli profiles 
of w and 9 do not change. On the other hand, \c(z)\ decreases (increases) near the plates 
(in the center of the fluid layer) so that the vertical profile of |c| flattens near the plates 
and becomes more peaked in the bulk near z — 1/2. The phase-lag of (p c (z) increases in the 
boundary layers with the through-flow. 

The above described flow-induced structural changes of c(z) can easily be understood 
within the L = approximation |32]j30[| to the concentration field balance (|3.1c|) 



Z{z)xi ^n _^L_ (4 , 8) 

in a critical TW that is propagating downstream with the critical frequency uf?. Thus for 
L = 0, i.e. in the absence of diffusive concentration currents — which, by the way, is quite 
a good approximation to the real situation of ethanol-water mixtures with L = 0.01 — the 
c(z)-profile is advectively slaved to the vertical velocity w(z) and the through- flow velocity 
U(z) = a P(z) Re: The prefactor i in fl4.8|) reflects the overall phase shift of 7r/2 between w 
and c. The effective z-dependent lateral velocity 

«?,/(*) = «?- a P(z)Jte (4.9) 
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that enters in the denominator of (|4.8|) renormalizes the profile w in the numerator of ( |4.8| ) 
so that the L = concentration profile 

cW!,I *lF^ (4l0) 

is determined by the quotient of these two velocity profiles. Since the phase velocity of the 
D-w&ve, = /k®, can well be approximated by 

Vc ~^c 0) +°$lpiRe (4. II) 

according to ( |4.3| ) with dtpi ~ 1.34 we can rewrite the lateral velocity that is effective for 
the concentration distribution profile ( 4.10|) of the D-wave as 



uf ff (z) ^^ 0) +aRe [1.34-62(1-2)]. (4.12) 

Thus, in the bulk u^Jz = 1/2) ps i>(°' — O.WaRe decreases with increasing through-flow 
while u®ff increases near the plates with Re, e.g., u®ff(z = 1/4) fa + 0.21aRe. This 
explains the flow-induced changes in the profile c (|4.10|) to be seen in Fig. |l^ for ip = —0.25 
with = 3.6. 

3. "Upstream" patterns 

In the absence of through-flow the field structure of TWs propagating to the left, i.e., 
in upstream direction is the mirror image of the Re = downstream patterns shown in the 
top part of Fig. |^. Therefore we have not included the Re = reference upstream pattern 
in Fig. |ll] where we show upstream patterns for Re = 0.3, 0.45, and 0.6. 

Switching on the through-flow decreases and eventually reverts the original Re = phase 
bending of the w and 6 waves (cf. right column in Fig. O) while the moduli \w\ and \9\ 



remain practically unaffected by the through-flow of Figs. [11] and [12[ On the other hand, 
the upstream concentration wave is significantly changed with increasing Re: \c(z) \ develops 
two side maxima in the upper and lower half of the fluid layer while flattening in the center 
part (cf. left column of Fig. [12] or also Fig. [Tl]). These structural changes of the c-wave 
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always occur in the .Re-range beyond the zero crossing of (Fig. |2|) where also and 
show a significant variation. Thus, we infer that these phenomena are related to each other. 

The flow-induced structural change of c(z) can be understood within the L = approx- 
imation 

(4.13) 



c(z) xnij—j- 

k c U eff( Z ) 



to the concentration balance (|3.1c|) in an upstream TW with effective lateral velocity 



«&/(*) 



v 



u 



aP(z) Re 



Again vf? = u^/k^ can well be approximated by 



-\vi 0) \ +a$tpx Re 



so that 



\vf ] \+aRe [1.34- 6 *(1 -z)\. 



(4.14) 



(4.15) 



(4.16) 



Note, however, that ignoring the dissipative contribution — LV 2 c on the Ihs of the concen- 
tration balance ( |3.1cj ) causes in the L = concentration profile ( [4.13| - [4.16] ) a divergence 
at 



2± = l±l 

2 2 



\ 



3 



via the zeros of u^jj(z) whenever 



Uo)| 
#e> 1 c 



1.34a 



(4.17a) 



(4.17b) 



The diffusive contribution in ( |3.1c|) prevents the divergence at z± leading to side maxima in 
\c(z)\ instead (cf. Fig. for Re = 0.45 and 0.6). Their location is shifted slightly towards 
the bulk in comparison with, e.g., Z-(Re = 0.45) = 0.1 and Z-(Re = 0.6) = 0.144. 

It should be noted that for the large Soret coupling ip = —0.25 the upstream fields 



diplayed in Fig. [XT], [12] for Re = 0.6 have not yet reached their asymptotic large- .Re-form - 
cf. also the bifurcation thresholds r^ 7 in Fig. |3[. For large Re the modulus \c(z)\ decreases 
substantially and becomes more flat, as was observed j|3| explicitly for smaller Soret coupling 
^ = -0.001, -0.01. 
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4- "Stationary" patterns 

Here we discuss the effect of through-flow on stationary patterns at a moderately positive 
Soret coupling, ip = 0.01, as a representative example. In the absence of through-flow, 



Re = 0, the perturbation fields (top part of Fig. 13 and full lines of Fig. FJ) are real, in 
phase, and laterally mirror symmetric around the vertical lines of maximal upflow (x = 0) 
and downflow (x = ±1). Thus, for Re = the lateral locations of largest alcohol surplus 
(deficiency) coincide with the largest vertical upflow (downflow) velocity that feeds the bulk 
with warm, alcohol rich (cold, alcohol poor) fluid from the bottom (top) plate. Note that 
for ip > the Soret effect causes concentration surplus (deficiency) at the warm (cold) plate. 
The critical convective concentration amplitude \c(z)\ is for Re = so large that the full 



line representing it in Fig. [14] lies outside the chosen plot range. The critical wavelength 
X^'(il> = 0.01) = 2.9833 of the Re = pattern (top part in Fig. [13]) is substantially larger 
than 2. 

Already a small through-flow changes the above described field structure of stationary 
convective perturbations dramatically. The wavelength decreases toward 2. The flow am- 
plitudes |u>(<2)| and even more conspicuously the concentration amplitude \c(z)\ decrease. 
The phase (p c (z) of the concentration field exhibits a strong vertical variation reflecting the 
almost passive advection by the flow that is roughly characterized by the parabolic lateral 
through-flow profile superimposed upon a roll-like closed-flow pattern. 

For ip < 0, however, amplitudes and phases of S"-patterns become more and more similar 
to those of -D-patterns when the two bifurcation branches r~ and r® approach each other 
with increasing Re in Fig. ||. However, the critical wave numbers of these two structures 
differ: For tfj = -0.001 we find £f = 1.001 and k s c = 0.998 and for ip = -0.01 the difference 
of the critical wave numbers (fcf = 1.089, k s c = 0.9678) increases. Thus, — kf increases 
with |^| when the bifurcation branches collide as described in Sees. [IVC4| and [IV C 5 . 
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V. ABSOLUTE AND CONVECTIVE INSTABILITY 

Whenever at a stability threshold the frequency is nonzero with a finite group velocity 

du{k) 



Vg " dk 



(5.1) 

k c 

one has to distinguish between spatiotemporal growth behavior of spatially extended and of 
spatially localized perturbations. The former having a form ~ e lkx have a positive growth 
rate above the bifurcation thresholds r c determined in the previous section. 

A. Wave packets, front propagation, and saddle point analysis 

A spatially localized perturbation, i.e., a wave packet superposition of plane wave ex- 
tended perturbations of a particular kind moves in the so called convectively unstable pa- 
rameter regime |P]-|7j with the velocity v g faster away than it grows — while growing in 



the frame comoving with v g the packet moves out of the system so that the basic conductive 
state is restored. In other words, the two fronts that join the wave packet's intensity enve- 
lope to the structureless state propagate both in the direction in which the packet center 
moves. On the other hand, in the so called absolutely unstable parameter regime the growth 
rate of the packet is so large that one front propagates opposite to the center motion. Thus 
the packet expands in the laboratory frame into the direction of packet motion as well as 
opposite to it [p|,p] . 

We should like to emphasize that we are dealing here only with a linear analysis of 
the convective fields. Thus we do not address the question whether, e.g., at a subcritical 
bifurcation the above decribed wave packet grows to a stable nonlinear state that then 
expands back into the system with a larger nonlinear front velocity so that the nonlinear 
structured state ultimately invades the region occupied by the homogeneous state. 

The boundary in parameter space between convective and absolute instability is marked 
by parameter combinations for which one of the fronts of the linear wave packet reverts 
its propagation direction in the laboratory frame: In the convectively unstable regime this 
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front propagates in the same direction as the center of the packet, in the absolutely unstable 
regime it moves opposite to it, and right on the boundary between the two regimes the front 
is stationary in the laboratory frame. This parameter combination can be determined by a 
saddle point analysis of the linear complex dispersion relation s(k) over the complex k plane 
p6| . Here we do not display the dependence of s on the control parameters r, Re, and ip. 

The condition of vanishing front propagation velocity is equivalent to finding the param- 
eters for which 

Us{k) = (5.2) 
with k denoting the appropriate saddle position of s(k) determined by solving 

«£>=0 (5.3) 
in the complex k plane (H6f . In Appendix [A 2| we describe our numerical method of finding 



the solution of ( |5.2| , |5.3| ). It yields the sought after surface in r — Re — ip parameter space, 
e.g., in the form of a function r c _ a (Re,ip) depending on Re and ip. This border r c _ a lies 
above the bifurcation threshold r c for growth of extended states. Thus for r < r c the basic 
conductive state is stable, for r c < r < r c _ a it is convectively unstable, and for r c _ a < r it 
is absolutely unstable. 

B. Ginzburg Landau amplitude equation approximation 



To solve (|5.2| , |5.3|) for r c _ a one has to determine the dispersion relation s(k; r, Re, ip) for 



complex k, i.e., one has to solve the eigenvalue problem Q3.9|) for complex k. This is in general 
a quite involved numerical task. A somewhat simpler, yet approximate method, is to use an 
expansion of s(k; r, Re, ip) that corresponds to approximate the full field equations by the 
Ginzburg Landau amplitude equation (GLE). We recapitulate the derivation of the relevant 
equations here for completeness. The comparison with results from the full field equation 
(Sec. |V G[ ) shows that the GLE yields quite useful approximations to the borderlines between 
absolute and convective instability. 
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1. Expansion of s(k;r) 

Consider Re and ip to be fixed for the moment so that we do not have to display them 
explicitly in the argument list of s. Under the assumptions (i) that the sought after saddle 
k in the complex k plane lies close to the critical wave number k c and (ii) that the relative 
distance 

Mc-a = — - 1 (5-4) 
r c 

between the convective/absolute border r c _ a and the critical Rayleigh number r c is small we 
expand 

+fi y r ^J + h - °- t - ■ ( 5 - 5 ) 
Here we have introduced for convenience the relative distance 

u, = — -1 (5.6) 

r c 

of the Rayleigh number r from its critical value r c for onset of convection. The higher order 
terms in ( |5.5|) should be of order /z 3 / 2 since for small < /i Cl only extended perturbations 
with (real) wave numbers out of a band of width k — k c ~ JJL can grow. 

2. Relation to linear amplitude equation 



The expansion coefficients of (|5.5| ) appear also in the linear parts of the complex GLE 



r (d t + v g d x ) A = n (1 + zc ) + il (1 + ic x ) d 2 



A 



+ nonlinear terms. (5-7) 
Here A(x,t) is the common complex amplitude of convection fields $ = (w,9,c) 

${x,z,t) = A{x,t) $(z) e i(kcX - wt) + c.c. (5.8) 
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that bifurcate out of the conductive state at /x = 0. The approximation ( |5.8| ) can be expected 
to be a good one as long as A is small and, more importantly, as long as the spatial field 
structure is well represented by that of the critical eigenf unctions 3>(z) e lkcX . 

The relations between the expansion coefficients of s and the coefficients in the amplitude 
equation are 



' ds\ . ( du " 

dk) = ~ l [ dk 



-iv g (5.9a) 



c 



\dk 2 J c r 
In ( |5.9cj ) we have used the relation [40 



' r ?*\ = r (?l-i^) - 1 + ic ° (5.9b) 

v dr J c \dr dr J r 

( Fp- ,\ of 2 

' > - 5o 'l + i Cl ). (5.9c) 



d 7 \ ( ^7 d r stab 



(5.10) 



\ dk 2 I \ dr dk 2 I 

\ / c \ / c 

to relate the second k derivative of the growth rate to the critical curvature £q = \ ( 9 ^2 ab ) 
of the marginal stability curve r sta b{k)- Note that the approximated dispersion relation ( |5.5| 
is precisely the one of the GLE approximation. 



3. Saddle and convective/absolute instability border 



With the notation ( |5.9f ) the saddle k ( |5T3"| ) of the approximated dispersion ( |5.5| ) lies at 



(§*) 

\9kJc 
(d 2 s\ 



k c 



l VnT 



1 + id 2£ 2 



(5.11) 



Then the condition (K3 



= 3? 



/i r 



9s' 
<9r 



1 

1 \dk 



( d 2 s\ 
\9k 2 J ( 



(5.12) 



yields the GLE approximation 



He 



?; 2 r 2 



(5.13a) 
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or 



r c _ a = (1 + fic-a) r c (5.13b) 
for the boundary between the convectively and absolutely unstable parameter regime 



4§-g8 



C. Convective and absolute instability against S, U, and D perturbations 



Here we present in Fig. [15] our numerical results for the borderlines between absolute and 
convective instability obtained from the full field equations (cf. Appendix |A2|) in comparison 
with the GLE results. 

There are three different types of extended perturbations in our systems, namely S, U, 
and D, against which the basic conductive state becomes unstable at the bifurcation thresh- 
old rf , , and that have been determined in Sec. |TV|. Consequently one has to investigate 
the spatiotemporal growth behavior of three different packets consisting of superpositions 
of S, U, or D plane wave perturbations. Their analysis along the lines of Sec. |V A| and 
Appendix |A 2| yields three functions r^_ a (circles), r^_ a (upwards pointing triangles), and 
r®_ a (downwards pointing triangles) each depending on Re and ip that mark the boundary 
surfaces in r — Re — ip parameter space between the convectively and absolutely regime of 
the basic state against type S, U, or D perturbations. 

In order to determine the three boundaries within the GLE approximation we have de- 
termined the derivatives (|5.9|) of the respective eigenvalues s s , s u , and s D at their respective 
critical points r c (Re,ip), k c (Re,ip) |43| . Then ( |5.13|) yields the functions r c _ a (i?e, tj)) for the 



three patterns. They are shown in Fig. [15| for a few ip as functions of Re by thick lines (S: 
dotted, U: full, D: dashed) together with the corresponding bifurcation thresholds r c (thin 
lines). For positive Soret coupling ip = 0.01, 0.001 one obtains within the GLE approxima- 
tion a local maximum of rf_ a (-0 = 0.01 at Re ~ 0.25, ip = 0.001 at Re ~ 0.1), when the 
product of r decreasing with Re and v g increasing with Re reaches a maximum. In the 
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ip = 0.01 case a second maximum is visible at small Re that is possibly associated with a root 
of c\. In contrast rf_ a of the full field equations (circles) is monotonously increasing with 
Re. For negative Soret coupling -0 = —0.001, —0.01, —0.1, and —0.25 r^_ a first decreases 
with increasing Re, coincides at its local minimum with r^ 7 when v g = 0, and afterwards 
increases with Re. The GLE approximation yields good quantitative agreement with the 
results of the full field equations for r^_ a (upwards pointing triangles) in the vicinity of the 
local minimum, while for higher Re the GLE results increase more strongly. For smaller 
Reynolds numbers this validity range is typically enlarged down to Re = for the ^-values 
presented here. 

As discussed in Sec. [IV C 5| we have limited our investigation of and r~ as functions 



of Re to cases, where the wave number gap in the D — S marginal stability curves has 
not yet appeared. Therefore, within the GLE approximation r®_ a and rf_ a are determined 
by the expansion coefficients of s ( |5.9|) at the critical Rayleigh number r c only as long as 
the calculation of rf? and rf was numerically possible. Nevertheless, for higher Re critical 
Rayleigh numbers still exist to the left and right hand side of the wave number gap (cf. 
Fig.|6]). We also found that the S saddle point of the full field equations (evaluated by 
succesively increasing Re) evolves monotonously when Re increases beyond the threshold 
where the wave number gap occurs. The wave number of this saddle point is always above 
k — 1 and increases with Re. 

The GLE results for r®_ a are increasing stronger than those of the full field equations 
(downwards pointing triangles). For ip = —0.001 the r®_ a stability limit seems to terminate 
close to the Reynolds number where the wave number gap (cf. Sec. |IV G 5| ) opens up in the 
D — S marginal stability curves. There seems to change sign when increasing Re further. 
A graphical analysis of this saddle that has moved into the lower complex fc-plane suggests 
that there 3fts < above a certain .Re-limit for all r so that this saddle can be ignored. 
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VI. CONCLUSION 

We have investigated the influence of an externally imposed horizontal shear flow on 
linear convective structure formation in binary fluid layers heated from below. To that 
end we have solved the linearized field equations for convective perturbations of the basic 
conductive state numerically with a shooting method. In addition we have checked our 
results — in particular some of the peculiar changes wih Re — by a variational calculation 
that gave a good agreement. We have determined for positive and negative Soret coupling ip 
the .Re-dependence of the critical bifurcation properties: stability thresholds, wave numbers, 
frequencies, and eigenfunctions for three different types of perturbations. The latter are 
identified by different characteristic exponents that cause perturbations to be stationary 
(S), downstream traveling (D), or upstream traveling (U) at Re = 0. 

The Hopf symmetry degeneracy of U and D perturbations at Re = and ip < is broken 
by a finite through- flow — wave numbers k^ and k^, frequencies and ou®, and bifurcation 
thresholds and develop differences. At small Re upstream traveling wave convection 
grows first upon increasing the heating since is depressed and r,? is shifted upwards. But 
then, with increasing through-flow the bifurcation lines and intersect giving rise to 
bistable bifurcation behavior. Eventually r]! flattens out and approaches for large Re the 
pure fluid's stability boundary rf (-Re, ip = 0) — a sufficiently strong shear flow eliminates the 
Soret induced coupling effects between the convective concentration field and temperature 
and velocity fields. This also holds for mixtures with positive ip. The bifurcation thresholds 
rf? monotonously curve upwards when increasing Re and collide with the rf threshold lines 
that for ij} < sharply drop downwards with growing Re. This behavior is easily understood 
by analyzing the .Re-variation of the paths of the relevant three eigenvalues in the complex 
7 — ui plane and of the variation of j(k, r). 

At ip > the stationary threshold rf rapidly approaches the tp — asymptote of a 
pure fluid from below — again the shear flow suppresses the vertical convective transport of 
Soret induced concentration gradients. Thus, any Soret effect's influence on the bifurcation 
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thresholds is eliminated by sufficiently large through-flow. 

Finally we have evaluated the borders r£_ oJ r®_ a , and rf_ a between convective and abso- 
lute instability for U, D, and S perturbations, respectively. To that end we have determined 
the relevant saddles of the dispersion relations s u (k), s D (k), and s s (k) in the complex wave 
number plane. These numerically exact results were compared with GLE approximations 
which agree reasonably well with the former. The latter were obtained by an expansion 
around the respective critical values and by evaluating all the coefficients that enter into the 
linear GLE for U, D, and S patterns. 

ACKNOWLEDGMENTS 

Support by the Deutsche Forschungsgemeinschaft is gratefully acknowledged. 



APPENDIX A: NUMERICAL ANALYSIS 

1. Shooting algorithm 

Here we describe our version of the shooting algorithm modified such as to better cope 
with the numerical problems caused by the concentration boundary layers near the horizontal 
plates. 

The equations fl3.1|) are written as an 8-D system of first order differential equations 

d z y - A{z, A) y = at 0<z<l . (Al) 



where the coefficient matrix follows directly from (|3.9|) . The eigenvalues 



A = [r stab (k), uj(k)] 7=0 or A = [ 7 (k), w(*)] r<tod (A2) 

are determined as a function of k either for 7 = on the marginal stability boundary r sta b{k) 
or for fixed r. The eigenf unctions 

y = (w, d z w, d 2 z w, dlw, 9, d z 9, C, 3,0 (A3) 
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fulfill the boundary conditions 



2/i = 2/2 = 2/5 = 2/s = at z = 0,l . 



(A4) 



following from (|3.3|). 

First (|A1[) is integrated with a fourth-order Runge-Kutta method with the initial vectors 



y 1 = e 3 , y 2 = e 4 , y 3 = e 6 , and y 4 = e 7 (A5) 

where is the unit vector in j-direction. Then an arbitrary linear combination 

y = E a, y> (A6) 

fulfills the boundary conditions at z = 0, but not at z = 1. The latter is guaranteed by 
choosing the correct coefficients aj. This leads to a set of homogeneous equations 



^ 1 2 3 4 ^ 

2/1 2/i 2/i 2/i 

1 2 3 4 

2/2 2/2 2/2 2/2 

12 3 4 

2/5 3/5 2/5 2/5 

y 2/8 2/1 2/1 2/s y 



«=i 



/ \ 

ai 

«2 

v ° 4 y 







(A7) 



The nontrivial solution of (|A7|) requires the determinant of the coefficients to vanish. Thus, 
this solvability condition but also the explicit evaluation of the eigenvalues requires the 
knowledge of the eigenf unctions at z — 1. This is a problem because of the strong variations 
in the boundary layers — see, e.g., Fig. |[ Our integrator is not able to shoot exactly 
from z = to z = 1 so that we can not evaluate the eigenvalues A with a fixed accuracy. 
Especially, this fact makes it difficult to find a good critical wave number. To understand 
this, we can expand the stability curve r sta b around k e c xact 



Ar = r stab (K xact + q)- r stab (K xact ) 



2 q 



d 2 r stab 



8k 2 



+ 0(q 3 ) 



(A8) 



and we find the defect 
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q 



IAr 



(A9) 



With £q ~ 0(1) the accuracy Ar strongly affects the accuracy of the critical wave number 
and of the coefficients entering the GLE. 

One could change the integrator. But the problem can be solved more easily. Using the 
mirror symmetry of the eigenfunctions (w, 8, () one only needs to shoot from z = to mid 
height z = 1/2. There the new boundary conditions are 

1 



2/2 = 2/4 = 2/6 = 2/8 = at z 



(A10) 



The new homogeneous equations replacing 



1 2 3 4^ 

2/2 2/2 2/2 2/2 

12 3 4 

yi vl 2/4 vl 

Vl Vl 2/1 2/e 

y vl 2/1 2/1 2/1 J 



are then 

/ \ 

ai 

a 2 
a 3 
y a 4 y 







(All) 



Now, shooting only from 2; = to z — 1/2 we can compute the eigenvalues with a fixed 
accuracy in our parameter range. A welcome secondary effect is that one saves CPU-time. 
However, it must be said that this trick only delays the boundary-layer-problem without 
eliminating it completely: If the parameters are such that the phase velocity \u c /k c \ reaches 
about 20 the boundary- layer at z = has become so thin, that it is no more possible to use 
our integrator. 



2. Saddle point analysis 

The borderline r c _ a between convective and absolute instability is defined by the condi- 
tions ( |5.2| , |5.3| ) on the dispersion relation after a continuation into the complex wave number 
plane, which also gives the saddle point position k and the frequency u(k). We have de- 
termined 3ls(k) and Qs(k) with a shooting method from z = to mid height z = 1/2 (see 
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Appendix [A 1|) for given r, Re, and if) and complex wave number fc. Then we solved for fixed 
Re and ip the nonlinear system of the three equations 

9te(«,r)=0, (A12) 



<99?s(k, r) 



(A13) 



(A14) 



that correspond to (^2,^3) with a Newton-Raphson method with backtracking. The solu- 



tion yields the borderline value r c _ a and the saddle position k of s in the complex fc-plane. 

The Jacobian matrix as well as the partial derivatives were obtained in discretized form 
by using central differences in the variables dtk, $sk, and r. Mainly two problems occur: one 
needs a good initial guess to ensure convergence of the Newton-Raphson method, and for 
higher Reynolds numbers, which are not discussed here, the shooting method starts to fail 
for the given accuracy limit (see Appendix A 1 ). 



We also used an alternative iterative method to evaluate the saddle points. It yields 

the same results as our first method, but it does not require partial derivatives. First we 
calculate the eigenvalue 

["(«), r(«)] 7=0 (A15) 

at a suitably chosen initial value k for the saddle position. Then, along a circle in the 
complex fc-plane of radius p around k we determine, say, n ~ 10 values 

["(A^.T^U-fl > J = 0,..,n-1 (A16) 

with 

% = S + 4c 0s ^V -infill • (A17) 



n J \ n 



The function values {u(n),uj(ka), .., u;(/c n _x)} are then used for a biquadratical fit 
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u(k) = ao + ai + a2 5i/c + a3 (3/c) 2 

+a 4 (^fc) 2 + a 5 SA;KA; (A18) 



and the analytically determined saddle of ( A18 ) is used as a new initial value k. Now, we 



reduce the circle radius p by a factor a and continue as above. We repeat this procedure 
until the changes of k have fallen below a suitable limit. This method is fast, robust, and 
easily programmable. However, the first p and the reducing factor a have to be chosen 
carefully. Typically we put p ~ 0(1CT 2 ) and a = y/2. 

APPENDIX B: COMPARISON WITH A VARIATIONAL METHOD 

In view of the unexpected Re dependence, say, of Fig. || we wanted to compare our results 
obtained with a shooting method with a completely different and independent method. To 
that end we used the variational method of Prigogine and Glandsdorff |]I9"| , |50|| . While this 
method does not yield exact results it still is very valuable in providing an independent 
check for our numerical analysis. 

The aim of this variational method is to find an functional, called local potential, and to 
minimize it. We start with the linearized field-equations 

d t u = — cr Re P d x u — d x p + aV 2 u 

-aRewd z P (Bla) 
d t w = —a Re P d x w — d z p + aV 2 w 

+a[(l+^)# + C] (Bib) 

d t 6 = -a RePd x 6 + V 2 6 + Raw (Blc) 

d t C = -aRePd x C + LV 2 C-tpV 2 e (Bid) 

= d x u + d z w . (Ble) 

Here u, w, 8, ( and p are convective disturbances from the conductive profiles and P(z) the 
Poiseuille shear flow ( 2.8b ). Following the procedure explained in Ref. we multiply the 
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above equations by increments of the respective fields: (|Bla| ) by —Su, ( |Blb|) by — Sw, (|Blc| ) 
by — 59 and (|Bld|) by — S(. Here, e.g., 



w = w° + Sw 



(B2) 



with w° being a solution of (|Bll ) that is not varied in the following. Then we add and obtain 

-\dt [(Su) 2 + (Sw) 2 + (S9) 2 + (SQ 2 
= d t u° Su + d t w° Sw + d t 6° S6 + d t (° S( 
+iT Re P d x u Su + a Re w d z P Su 
+d x p Su — a Su V 2 u + a ReP d x w Sw 
-a Sw [(1+ %l))9 + C] - ° SwV 2 w - RawS9 
-SO V 2 6 + aReP dj SO + d z p Sw 

-L5(V 2 ( + ip8(V 2 e + aRePd x (8( (B3) 



when making use of relations like 



1 



— d t w5w = —d t (w° + Sw) Sw = —-d t (Sw) 2 — d t w° Sw 



(B4) 



We expand the fields in plane waves, $ = $(z)e e st , and we integrate over the entire x — z 



crosssection of the fluid layer. Following closely the prescriptions of Ref. [50] we then otain 
the local potential 



dz 



k 2 



d z w°d z w+s (w°w + e°e + c°c 



+a d z w° d z w 



a 



, } , {dlw) 2 + ^d 2 z w^lw 



+ikaReP(w°w + 6° 6 + C°C) + - (d, 



,w) 



+ak 2 (w°w + -e°e + - 

-Ra w*6+ l - (d z 6) 2 + - (d z () 2 -a(°w 



[l + ^)6°w 



-V> k 2 9° ( + ij d 2 z 9° ( - r a Reid d z w d z P 

k 



+ - a ReP d z w° d z w 
k 



(B5) 
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Here we have elimitated the pressure by using ( |Bla| ) 

k z k 

- l -o ReuP d z P - ^dlw° + ad z w° . (B6) 

Furthermore, partial integrations in z have been applied. One can check directly that the 
Euler-Langrange equation for, say, ( 

S * S ^ -0 (B7) 



d( dz g h£ 
leads to 

[s + ikaReP] (° = L (d 2 ( - k 2 (°) - V (d 2 -k 2 )P . (B8) 
This is together with the a-posteriori subsidiary condition Q = C° the equation ( [B 1 d|) . An 



analogous calculation leads with respective subsidiary conditions 6 = 9°, w = id to the 
other differential equations. 

To get the critical values we expand the unknows w, 9 and ( 

Ni Ni 

= E ^ fi (*) ®° W = E (B9a) 

3=1 3=1 

= S°(*) = E ( B9b ) 

3=1 3=1 

?(*) = E 5 hj(z) (°(z) =J2(j hi{z) (B9c) 

3=1 3=1 

in functions /j, gj, and /ij that satisfy the NSI boundary conditions. As trial functions we 
use 

f s (z) = [z(l - z)) 2j (BlOa) 

9j (z) = z(l -z)(2z- 1) 2{ ' J - 1] (BlOb) 

hj(z) = [z(l - z)] 2 ^ . (BlOc) 



We insert ( |B"9"| , |B1(J| ) into the expression for we minimize with respect to the variational 



parameters Wj,9j and Q, and we use the subsidiary conditions 
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w 3 =w^ 63 = 6% = C° • (Bll) 

Then we obtain a system of N = N\ + A2 + A3 linear homogeneous equations Ax. = where 
the vector x has iV components containing w®, 9®, and (®. For nontrivial solutions to exist 
the determinant of A must be zero, so that we have to solve the equation 

det A[r stab (k), oo(k)] J=0 = . (B12) 

With Ni = 1, N 2 = 3, iV 3 = 4 the software package MATHEMATICA is able to solve ( gflg ) 
exactly. 

The variational results display the same structural properties as the shooting results and 
they agree within less than 1% — 2% for Re 7^ 0. For Re = the agreement is better. We 
found, e.g., for tjj = -0.25 the critical values r c = 1.33497, k c = 1.00451, and uj c = ±11.20027 
differing by less than 0.2% from the shooting results. As an interesting aside we mention 
that the variational calculus yields in the absence of through-flow ip^ = —Lj (1 + L) exactly. 
Furthermore, the separation ratio ipo for which k c = is obtained as ipo = L/(f — L) with 
/ = 6695603/25739142 = 0.2601... which has to be compared with the analytically exact 
result / = 34/131 = 0.2595... 0,0. 
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FIGURES 

FIG. 1. Eigenvalues for small Reynolds numbers. The variation of growth rates 7 and frequen- 
cies ijj with Re resulting from an expansion up to linear order in Re is shown schematically. The 
lifting of the symmetry degeneracy, = s[ , by the through-flow and the behavior of S3 is 



described in Sec. Ill D 



FIG. 2. Evolution of Hopf bifurcation properties with through-flow. Shown are critical Rayleigh 
numbers (a), reduced wave numbers (b), and frequencies (c) for upstream (full lines) and down- 
stream (dashed lines) traveling perturbations that are symmetry degenerate for Re = 0. Parameters 
are L = 0.01, a = 10, ^ = -0.1. 

FIG. 3. -Re-dependence of bifurcation thresholds for negative Soret coupling ip. The stability 
boundaries rf (dotted lines), rj? (dashed lines), and (full lines) are shown for some representative 
■0-values as indicated. The behavior beyond the Re values where r£ and r~ collide is discussed in 
Sec. IV C 5| . Parameters are L = 0.01, a = 10. 



FIG. 4. (a) Frequencies tu{r) and (b) growth rates j(r) of the two eigenvalues that cause 
collision of the stability curves rj? and in Fig. |3[ Thick (thin) lines and symbols referring to 
Re = {Re = 0.4) were obtained for a fixed wave number k = 1. Arrows indicate deformation 
directions of the curves. The Re- variation of the zeros of 7 at r u (upwards pointing triangle), 
r D (downwards pointing triangle), and r s (circle) is discussed in Sec. [IV C 4 . Parameters are 
L = 0.01, a = 10, tp = -0.01. 

FIG. 5. Schematic variation of the eigenvalues of Fig. || in the complex 7 — w-plane. Arrows 
indicate the motion of the eigenvalues with increasing r. Line styles and symbols are those of 

Fig.a 



February 9, 2008 48 

FIG. 6. The "mountain landscape" of ^ D l s over the k—r plane for Re = 0.4 (a) and Re = 0.475 
(b). Gray-scales show the height of ^ D l s in the range where 7 > 0. Thick lines are marginal 
stability curves where 7 = 0. In the white parts of the figures ^ D / s < 0. Dashed lines are iso-7 
lines. The zeros of j D / s marked by circle and downwards pointing triangle are those of Fig. |]b. 
Parameters are L = 0.01, a = 10, and i/i = —0.01. 

FIG. 7. (a) Critical Rayleigh numbers and (b) wave numbers for positive ip versus through-flow 
Reynolds numbers. Parameters are L = 0.01, a = 10. 

FIG. 8. To give an impression of the bifurcation surfaces in r — Re — ip space we show r\! (thin 
full lines), (thin dashed lines), and rf (thin dotted lines) together with their ^-dependence for 
Re = (thick lines). Parameters are L = 0.01, a = 10. 

FIG. 9. Spatial structure of D patterns that propagate downstream, i.e., to the right. Contours 
of the fields w, 9, c in vertical cross sections of the fluid layer are shown for different Re as 
indicated. White (black) implies large (small) field values. The largest vertical up (down) flow is 
at x = 0(rbl). Contour lines mark fractions ±(0, 1,2, 3,4) /5 of the field maxima. Parameters are 
L = 0.01, a = 10, tfi = —0.25. The critical frequencies are about w 11.21 + 41.9 Re. 

FIG. 10. Vertical profiles of D patterns. Shown are moduli (left column) and phases (right 
column) of the critical complex amplitudes that are symmetric around the mid position z = 1/2. 



Their normalization is explained in Sec. IV F. Parameters are L = 0.01, a = 10, ip = —0.25. 



FIG. 11. Spatial structure of U patterns. They propagate for Re — > upstream, i.e., to the left. 
Contours of the fields w, 0, c in vertical cross sections of the fluid layer are shown for different Re 
as indicated. White (black) implies large (small) field values. The largest vertical up (down) flow 
is at x = 0(±1). Contour lines mark fractions ±(0, 1,2, 3,4)/5 of the field maxima. Parameters 
are L = 0.01, a = 10, ip = —0.25. The critical frequencies are about uj]! as —11.21 + 41.9 Re. 
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FIG. 12. Vertical profiles of U patterns. Shown are moduli (left column) and phases (right 
column) of the critical complex amplitudes that are symmetric around the mid position z = 1/2. 



Their normalization is explained in Sec. IV F. Parameters are L = 0.01, a = 10, ip = —0.25. 



FIG. 13. Spatial structure of S patterns. They are stationary for Re — ► 0. Contours of 
the fields w, 0, c in vertical cross sections of the fluid layer are shown for different Re as indi- 
cated. White (black) implies large (small) field values. The largest vertical up (down) flow is at 
x = 0(±1). Contour lines mark fractions ±(0, 1, 2, 3, 4)/5 of the field maxima. Parameters are 
L = 0.01, a = 10, if) = 0.01. The critical frequencies are about ~ 41.9 Re. 

FIG. 14. Vertical profiles of S patterns. Shown are moduli (left column) and phases (right 
column) of the critical complex amplitudes that are symmetric around the mid position z = 1/2. 



Their normalization is explained in Sec. IV F. Parameters are L = 0.01, a = 10, if) = 0.01. 



FIG. 15. Borderlines r c _ a between absolute and convective instability versus through-flow rate. 
Symbols and thick curves represent r c _ a obtained from the full field equations and from the GLE, 



respectively. Thin curves show stability curves r c discussed in Sec. |V). The different types of 
perturbations are identified by circles and dotted lines (5), upwards pointing triangles and full 
lines (U), and downwards pointing triangles and dashed lines (D). The case of a pure fluid, 
ip = 0, is shown by dash-dotted lines: thick one for r^_ a and thin one for r~. Parameters are 
L = 0.01, a = 10, and if) as shown. 
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Figure 6 
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Figure 9 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 




Figure 10 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 



CD 







^^^^^ 


• 


1131 


i \ \ V 



Re=0.45 



CD 



CD 




Figure 11 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 




Figure 12 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 



Re=0.0 


-0 = 0.01: Stationary (S) 


» 




■ 

Re=0.15 


J 


• 




• 




T> -ii i: 






• 







-1 1 



Figure 13 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 




Figure 14 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 




Figure 15 

Influence of through-flow on linear pattern formation properties 
in binary mixture convection 

Ch. Jung et al., Phys.Rev. E 



